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Abstract. We consider a random system of equations Xi + Xj = 6(j j\(mod 2), 
(xu G {0,1}, &(u,i;) = ^(t;,u) £ {0,1}), with the pairs {i,j) from E, a symmetric 
subset of [n] x [n] . E is chosen uniformly at random among all such subsets of a given 
cardinality m; alternatively (i, j) G E with a given probability p, independently 
of all other pairs. Also, given E, Pr{6e = 0} = Pr{be = 1} for each e £ E, 
^^ ' independently of all other bg/ . It is well known that, as m passes through n/2 

r) , (p passes through 1/n, resp.), the underlying random graph G(n, #edges = m) , 

{G{n, Pr(edge) = p) , resp.) undergoes a rapid transition, from essentially a forest 
of many small trees to a graph with one large, multicyclic, component in a sea of small 
C^ , tree components. We should expect then that the solvability probability decreases 

precipitously in the vicinity of ?7i ~ n/2 (p ~ I/jt.), and indeed this probability is 
of order (1 — 2m/n)^''^ , for m < n/2 ((1 — pn)^''^ , for p < 1/n, resp.). We show 
that in a near-critical phase m = (n/2)(l + Xn~^'^) (p = (1 + Xn~^'^)/n, resp.), 
A = o(n^/^^), the system is solvable with probability asymptotic to c{X)n~^'^'^ , 
for some explicit function c(A) > . Mike MoUoy noticed that the Boolean system 
with 6e = 1 is solvable iff the underlying graph is 2 -colorable, and asked whether this 
connection might be used to determine an order of probability of 2-colorability in the 
near-critical case. We answer Mike's question affirmatively and show that probability 
of 2-colorability is < 2-l/4el/8c(A)n-l/l2 ^ and asymptotic to 2-l/4el/8c(A)n-l/l2 
at a critical phase A = 0(1) , and for A -^ —oo . (Submitted to Electronic Journal of 
Combinatorics on September 7, 2009.) 
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k> , 1. Introduction. A system of 2 -linear equations over GF{2) with n Boolean 

\-{ ' variables xi, . . . , x„ E {0, 1} is 



C^ 



(1.1) Xi + Xj = bij (mod 2), bij = bj^i £ {0, 1}; {i ^ j). 

Here the unordered pairs (i, j) correspond to the edge set of a given graph G on 
the vertex set [n] . The system (1.1) certainly has a solution when G is a tree. It 
can be obtained by picking an arbitrary Xi S {0, 1} at a root i and determining 
the other Xj recursively along the paths leading away from the root. There is, 
of course, a twin solution Xj = 1 — Xj, j G [n] . Suppose G is not a tree, i.e. 
i{G) := e{G) — v{G) > 0. If T is a tree spanning G , then each of additional edges 
ei, . . . , e^(G)_|_i forms, together with the edges of T, a single cycle Gt,t< i{G) + 1 . 
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2 BORIS PITTEL AND JI-A YEUM 

Obviously, a solution Xj{T) of a subsystem of (1.1) induced by the edges of T is a 
solution of (1.1) provided that 

(1-2) bij =Xi{T) + Xj{T), (i,j) = ei,...,e£(G)+i; 

equivalently 

(1.3) Y^ 6^ = (mod 2), t = l,...,e{G) + l. 

e&E(Ct) 

So, intuitively, the more edges G has the less likely it is that the system (1.1) has 
a solution. We will denote the number of solutions by S{G) . 

In this paper we consider solvability of a random system (1.1). Namely G is 
either the Bernoulli random graph G{n,p) = G{n, Pr(edge) = p) , or the Erdos- 
Renyi random graph G{n, m) = G{n, ^ of edges = m) . Further, conditioned on the 
edge set E{G{n,p)) {E{G{n,m) resp.), 6e 's are independent, and Pr(6e = 1) =p, 
for all e. We focus on p = 1/2 and p = 1. p = 1/2 is the case when ftg 's are 
"absolutely random". For p= 1, 6e 's are all ones. Mike Molloy [17], who brought 
this case to our attention, noticed that here (1.1) has a solution iff the underlying 
graph is bipartite, 2-colorable in other words. 

It is well known that, as m passes through n/2 {p passes through 1/n, resp.), the 
underlying random graph G{n,m), {G{n,p), resp.) undergoes a rapid transition, 
from essentially a forest of many small trees to a graph with one large, multicyclic, 
component in a sea of small tree components. Bollobas [4], [5] discovered that, for 
G{n, m) , the phase transition window is within [mi , 7712] , where 

mi,2 =n/2±An2/3, A = e(ln^/2n). 

Luczak [14] was able to show that the window is precisely [?tii,?t22] with A — ?• 00 
however slowly. (See Luczak et al [16], Pittel [19] for the distributional results 
on the critical graphs G{n,m) and G{n,p).) We should expect then that the 
solvability probability decreases precipitously for m close to n/2 {p close to 1/n 
resp.). Indeed, for a mw/tigraph version of G{n,m), Kolchin [13] proved that this 
probability is asymptotic to 

^^^ (l-(l-2p)7)i/4' ^- „' 

if limsup7 < 1. See Creignon and Daude [9] for a similar result. Using the results 
from Pittel [19], we show (see Appendix) that for the random graphs G{n,^n/2) 
and G{n,p = 7/n) , with lim sup 7 < 1 , the corresponding probability is asymptotic 
to 



(1-7) 



1/4 



''^'^ (l-(l-2ph)V.°» 



„2 



-p+—p{l-p) 



The relations (1.4), (1.5) make it plausible that, in the nearcritical phase |m— n/2| = 
0(n^''^), the solvability probability is of order n~^'^'^ . Our goal is to confirm, 
rigorously, this conjecture. 
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To formulate our main result, we need some notations. Let {fr}r>o be a sequence 
defined by an implicit recurrence 



(1.6) 



/o — 1) /_^ fkfr-k — ^r, ^r '■- 



(dry. 



fe=0 



25r32r(3^)!(2r)!' 



Equivalently, the formal series '^^x'^fr, "^^x^^r (divergent for all x ^ 0) satisfy 



:i.7) 



E-v. =E 






It is not difficult to show that 

e 



1. 



2 v'-9-^^-T' ^>0- 



For y, A G M, let A{y, A) denote the sum of a convergent series, 



(1.9) 



Aiy,X) 



-AV6 



E 



(i3^/3A)' 



3(i'+i)/3 ^ fc!r[(y + 1 - 2k) /3] ■ 



We will write B^ ~ C„ if lim„^oo Bn/C'n = 1, and B^ < Cn if limsup„ Bn/C'n < 
1. Let Sn denote the random number of solutions of (1.1) with the underlying 
graph being either G{n,m) or G{n,p), i. e. Sn = S{G{n,m)) or Sn = S{G{n,p)) , 
and the (conditional) probability of 6e = 1 for e € E{G{n,m)) (e G E{G{n,p)) 
resp.) being equal p. 

Theorem 1.1. (i) Let p = 1/2. Suppose that 

l + An-i/3 



n 



m = — (1 + An -"^Z^), p 



(1.10) 

Then, for both G{n,m) and G{n,p), 



o[n 



1/12> 



n 



(1.11) 
where 

(1.12) 



Pr(5„ > 0) ~ n-^/'^ciX), 



/> 



g3/8(2^)i/2 ^ ^ ^(o_25 + 3r, A), AG (-oo, oo); 



r>0 



c(A) := <^ 



=3/8|X|l/4 



=3/8 



L 4 • 33/4 



A^/^exp(-10AV81), 



A- 

A — > 00. 



00; 



(ii) Let p=l. Then, with c(A) replaced by ci(A) := 2-i/''e^/*c(A) , ("i.P; holds for 
both G{n,m) and G{n,p) if either A = 0(1), or A — t- —00, |A| = o{n^'^'^) . For 
A ^00, A = o(ni/i2), 

Pr(5. >0)<n-i/i2^i(A). 
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Notes. 1. For G{n,m) with A — t- — oo, and p = 1/2, our result blends, 
qualitatively, with the estimate (1.4) from [13] and [9] for a subcritical multigraph, 
and becomes the estimate (1.5) for the subcritical graphs G{n,m) and G{n,p). 

2. The part (ii) answers Molloy's question: the critical graph G(n,m) {G{n,p) 
resp.) is bichromatic (bipartite) with probability -^ ci{\)n~^'^'^ . 

Very interestingly, the largest bipartite subgraph of the critical G{n,p) can be 
found in expected time 0{n), see Coppersmith et al [8], Scott and Sorkin [21] and 
references therein. The case A — t- oo of (ii) strongly suggests that the supercritical 
graph G{n,p = c/n), {G{n,m = cn/2) resp.), i. e. with liminfc > 1, is bichro- 
matic with exponentially small probability. In [8] this exponential smallness was 
established for the conditional probability, given that the random graph has a giant 
component. 

Here is a technical reason why, for A = 0(1) at least, the asymptotic probability 
of 2 -colorability is the asymptotic solvability probability for (1.1) with p = 1/2 
times 2~'^'^e^'^ . Let Ci{x) (C|(x) resp.) denote the exponential generating func- 
tions of connected graphs G (graphs G without odd cycles resp.) with excess 
e{G) — v{G) = £ > 0. It turns out that, for |x| < e~^ (convergence radius of Ce{x) , 
C^{x)), and x — )■ e~^. 



Gf(x) 



^Ceix), i>0, 

JCo(x) + ln(2-i/4ei/8)+o(l), i = 0. 



Asymptotically, within the factor e "v ^ ' , this reduces the problem to that 
for p = 1/2. Based on (1.5), we conjecture that generally, for p € (0, 1] , and the 
critical p, Pr(S'„ > 0) is that probability for p = 1/2 times 



:i.ll) (2^)-V4exp 



i^-pr , 1 



+ 7T 



(For ^ = 0, Pi{Sn > 0) = 1 obviously.) 

3. While working on this project, we became aware of a recent paper [10] by 
Daude and Ravelomanana. They studied a close but different case, when a system 
of m equations is chosen uniformly at random among all n(n — 1) equations of 
the form (1.1). In particular, it is possible to have pairs of clearly contradictory 
equations, Xi + Xj = and Xi + Xj = 1. For m = 0{n) the probability that none 
of these simplest contradictions occurs is bounded away from zero. So, intuitively, 
the system they studied is close to ours with G = G{n,m) and p = 1/2. Our 
asymptotic formula (1.9), with two first equations in (1.10), in this case is similar 
to Daude-Ravelomanana's main theorem, but there are some puzzling differences. 
The exponent series in their equation (2) is certainly misplaced; their claim does 
not contain our sequence {fr}- 

As far as we can judge by a proof outline in [10], our argument is quite different. 
Still like [10], our analysis is based on the generating functions of sparse graphs 
discovered, to a great extent, by Wright [23], [24]. We gratefully credit Daude 
and Ravelomanana for stressing importance of Wright's bounds for the generating 
function Gi{x) . These bounds play a substantial role in our argument as well. 
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4. We should mention a large body of work on a related, better known, 2 — SAT 
problem, see for instance BoUobas et al [6], and references therein. It is a problem 
of existence of a truth-satisfying assignment for the variables in the conjunction 
of m random disjunctive clauses of a form Xi V Xj, {i,j G [n]). It is well known, 
Chvatal and Reed [7], that the existence threshold is m/n = 1. It was proved in 
[6] that the phase transition window is [mi, 1712], with 

"^1,2 =n±An', |A|— T-oo however slowly, 

and that the solvability probability is bounded away from both and 1 iff ttt, = 
n + 0(n2/3). 

5. A natural extension of the system (1.1) is a system of fe -linear equations 
(1.12) J]x, = 6e(mod2), 

where e runs over a set E of (hyper)edges of a /;; -uniform hypergraph G , k > 2, on 
the vertex set [n] , Kolchin [13]. Suppose G is chosen uniformly at random among 
all A: -uniform graphs with a given number m of edges, and, given G , the beS are 
independent Bernoullis. It will be interesting to study, for k > 2, the limiting 
solvability probability as a function of m/n. See [13] for some thought-provoking 
results on the behavior of the number of hypercycles in this random hypergraph. 

The paper is organized as follows. In the section 2 we work on the G{n,p) and 
p = 1/2 case. 

Specifically in the (sub)section 2.1 we express the solvability probability, Pr(S'n > 
0) , and its truncated version, as a coefficient by x" in a power series based on the 
generating functions of the sparsely edges (connected) graphs. We also establish 
positive correlation between solvability and boundedness of a maximal "excess", 
and determine a proper truncation of the latter dependent upon the behavior of 
A. In the section (2.2) we provide a necessary information about the generating 
functions and their truncated versions involved in the formula and the bounds 
for Pr{Sn > 0) . In the section 2.3 we apply complex analysis techniques to the 
"coefficient by x" " formulas and obtain a sharp asymptotic estimate for Fr{Sn > 
0) for |A| =o(ni/i2). 

In the section 3 we transfer the results of the section 2 to the G{n, m.) and 
p = 1/2 case . 

In the section 4 we establish the counterparts of the results from the sections 2,3 
for G{n,p) , G{n, m) with p = 1. An enumerative ingredient of the argument is an 
analogue of Wright's formulas for the generating functions of the connected graphs 
without odd cycles. 

In Appendix we prove some auxilliary technical results, and an asymptotic for- 
mula for Pr{Sn > 0) in the subcritical case, i. e. when the average vertex degree 
is less than, and bounded away from 1 . 

2. Solvability probability: G{n,p) and p= 1/2. 

2.1. Representing bounds for Pr(S'„ > 0) as a coefficient of x" in a 
power series. 
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Our first step is to compute the probability of tlie event {Sn > 0} , conditioned 
on G{n,p). Given a graph G = {V{G),E{G)) , we denote v{G) = \V{G)\, e(G) = 
\EiG)\. 

Lemma 2.1.1. Given a graph G on [n], let c{G) denote the total number of its 
components Hi . Then 



c{G) e{H,)-{v{H,)-l) 

Pr(5„>0|G(n,p) = G) = J] - 



i = l 



X(G) 



, X{G):=e{G)-n + c{G). 



Consequently 



Pr{Sn > 0) = E 



X(G(n,p)y 



Proof of Lemma 2.1.1. Recall that, conditioned on G{n,p) , the edge variables 
be are mutually independent. So it is suffices to show that a system (1.1) for a 
connected graph H , with independent b^, e £ E{H) , such that Pr(6e = 1) = 1/2, 
is solvable with probability (1/2)^"'"-^ , where £ = e{H) — v{H) . 

Let T be a tree spanning H . Let x(T) := {xi{T)}i^y(^}j-^ be the solution of the 
subsystem of (1.1) corresponding to v{H) — 1 edges of T, with x^j, = 1 say, for a 
specified "root" i^. x(T) is a solution of the whole system (1.1) iff 



(2.1.1) 



Xi{T) + Xj{T), {{ij) = e), 



for each of e{H) - {v{H) - 1) = ^ + 1 edges e G E{H) \ E{T) . By independence 
of be 's, the probability that, conditioned on {&e}ee_E(T) ? the constraints (2.1.1) are 
met is (1/2)^+^. (It is crucial that Pr(6e = 0) = Pr(6e = 1) = 1/2.) Hence the 
unconditional solvability probability for the system (1.1) with the underlying graph 
H is (1/2)^+1 as well. D 

Note. For a cycle G C H, let be = X]ee_B(c)^e' '^^^ conditions (2.1.1) are 
equivalent to be being even for the i + 1 cycles, each formed by adding to T 
an edge in E{H) \ E{T). Adding the equations (1.1) over the edges of any cycle 
C (^ H , we see that necessarily be is even too. Thus our proof effectively shows 
that 

Pr ^ (I {6c is even} 

CCH 



Using Lemma 2.1.1, we express P{S{n,p) > 0) as the coefficient by x"" in a 
formal power series. To formulate the result, introduce Gi^x), the exponential 
generating function of a sequence {G{k, k + £)}k>i , where C{k, k + i) is the total 
number of connected graphs H on [k] with excess e{H) — v{H) = I. Of course, 
G{k,k-V€) = unless -\<i< (3) -k. 
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\T.[f,\C'i-) 



i>-i 



Lemma 2.1.2. 

(2.1.2) Pr{Sn > 0) =N{n,p) [x"] exp 

(2.1.3) N{n,p):=nlq-'/'(^ 



Proof of Lemma 2.1.2. The proof miniicks derivation of the "coefficient-of 
x'^- expression" for the largest component size distribution in [19]. 

Given a = {ak/} , such that '^^k i^^k/ = n, let Pn{oi) denote the probability 
that G{n,p) has ak,e components H with v{H) = k and e{H) — v{H) = i. To 
compute Pnioi) , we observe that there are 

n! 



n(A:!)°'=-^afe,^! 
k,e 



ways to partition [n] into ^j, ^ Ok/ subsets, with ak/ subsets of cardinality k and 
"type" i. For each such partition, there are 

Y[[c{k,k+i)r'-' 

k,e 

ways to build ak,e connected graphs H on the corresponding ak/ subsets, with 
v{H) = k, e{H) — v{H) = i. The probability that these graphs are induced 
subgraphs of G{n,p) is 



n[: 

k,e 



fe+^o(2)-(fe+^)' 



p q 



P 

73/2 



n 

k,e L 



P] ,^V2 



^^ Sfc e ^ ^k,e = n. The probability that no two vertices from two different subsets 
are joined by an edge in G{n,p) is q^ , where r is the total number of all such pairs, 
i. e. 

r=^k'^( I'j + n Yl kik2ak,,e,ak^,e2 



kl 



(fcl,£l)/(fc2,«2) 



k,£ \ k,e J 

1 v^ , 2 n"^ 

-^k ak,i + —. 



k,e 



Multiplying the pieces. 



Pr,ic)=N{n,p)l[ 



1 



ki "'^'^^ L 



{p/qyC{k,k + i) 

k\ 
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So, using Lemma 2.1.1, 

(2.1.4) Pr(5„>0)=iV(n,p)^[] — 



a kJ 



Oik,. 



[l/2f+\p/qYC{k,k + i) 
k\ 



ak,l 



Notice that dropping factors (1/2)^+-'^ on the right, we get 1 instead of Pr(5„ > 
0) on the left, i.e. 



(2.1.5) 



l = iV(n,p)J]n 



OL k,i 



Oik/- 



ip/qYcik,k + e) 
k\ 



So, multiplying both sides of (2.1.4) by ^f . and summing over n > 0, 
(2.1.6) 



E^ 



Pr(5n > 0) 
N{n,p) 



E n 

y^ ka.k^i<oo k,i 



^kctk 



Oik/- 



[l/2Y+\p/qYC{k,k + l) 
k\ 



-I CtkJ 



:exp 



:exp 



e k 



^^^ C{k,k + l)x^ 



k\ 



We hasten to add that the series on the right, whence the one on the left, converges 
for X = only. Indeed, using (2.1.5) instead of (2.1.4), 



(2.1.7) 



exp 



E(p/9)'C'(^) 



E^ 



1 



xq' 



3/2 



N{n,p) 






n\ q 



n2/2 



OO, 



for each x > 0. Therefore, setting p/2q = pi/qi , (^i = 1 — Pi ), 

Y,{p/'^qYc,{x) = Y,{pihiYc,{x) = OO, vx > o, 

as well. D 

Note. Setting pjq = w , x = yw , in (2.1.7), so that p = w/{w + 1), q = 
l/{w + 1) , we obtain a well known (exponential) identity, e. g. Janson et al [12], 



exp 



Y^ w'^Ceiyw) 



i>-i 



Y^l^{w + iY^) 

n>0 



the right expression (the left exponent resp.) is a bivariate generating function 
for graphs (connected graphs resp.) G enumerated by f(G) and e(G). Here is a 
similar identity involving generating functions of connected graphs G with a fixed 
positive excess. 



(2.1.e 



exp 



Y.w'G,{x) 



Y^w'^Er 



e>i 



r>0 
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where Eq{x) = 1, and, for £ > 1, Ei{x) is the exponential generating function 
of graphs G without tree components and unicycHc components, that have excess 
i{G) = e{G) — v{G) = l, see [12]. In the hght of Lemma 2.1.2, we wih need an 
expansion 



(2.1.9) exp 



\Y.w'C,{x) 



2 



S^wTrix). 



r>0 



Like Er{x) , each power series Fr{x) has nonnegative coefficients, and converges 
for |x| < e~^ . 

By Lemma 2.1.2 and (2.1.8), 

Pr(5„ > 0) =Nin,p)Y,(j-) [^"] {e^^"^^r(x)} ; 
(2.1.10) '■^o ^ ^^ 

H{x):=^C.i{x) + Igo{x). 
p 2 

Interchange of [x"^] and the summation is justifiable as each of the functions on 
the right has a power series expansion with only nonnegative coefficients. That is, 
divergence of X]^(p/2(7)^Q(x) in (2.1.6) does not impede evaluation of Pr(S'„ > 0) . 
Indirectly though this divergence does make it difficult, if possible at all, to obtain 
a sufficiently sharp estimate of the terms in the above sum for r going to oo with 
n, needed to derive an asymptotic formula for that probability. Thus we need to 
truncate, one way or another, the divergent series on the right in (2.1.6). One of the 
properties of Ge{x) discovered by Wright [23] is that each of these series converges 
(diverges) for \x\ < e~^ (for \x\ > e~^ resp.). So, picking L >0, and restricting 
summation range to £ € [—1, L] , we definitely get a series convergent for \x\ < e~^ . 
What is then a counterpart of Pr(S'„ > 0)? Perusing the proof of Lemma 2.1.2, 
we easily see the answer. 

Let G be a graph with components Hi,H2, .... Define £{G) , a maximum excess 
of G, by 

£{G) = max[e{H,)-viHi}]. 

i 

It can be easily seen that £{G) is monotone increasing, i. e. £{G') < £{G") if 
G' CG". Let £n = £{G{n,p)). 



Lemma 2.1.3. 

(2.1.11) Pr{Sn >0,£„<L)= N{n,p) [x"] exp 



^1.(1'""^' 



The proof of (2.1.11) is an obvious modification of that for (2.1.2). 

If, using (2.1.11), we are able to estimate Pr(S'„ > 0, <?„ < L) , then evidently we 
will get a lower bound of Pr(S'„ > 0) , via 

(2.1.12) Pr(5„ > 0) > Pr(5„ > 0, £n < L). 

Crucially, the events {Sn > 0} and {£n < L} are positively correlated. 



10 
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Lemma 2.1.4. 

(2.1.13) 



Pr{S^ > 0) < 



Pr{Sn > 0, gn < L) 
Pr{£n < L) ■ 



Note. The upshot of (2.1.12)-(2.1.13) is that 

Pr(5„>0)- Pr(5„>0,£:„<L), 

provided that L = L{n) is just large enough to guarantee that Pr(i£'„ < L) 
Proof of Lemma 2.1.4. By Lemma 2.1.1, 



Pr(5„ > 0, £:„ < L) = E 



XiG{n,p)) 



^{£(G(n,p))<L} 



where X{G) = e{G) — n + c{G) . Notice that (l/2)'''"(^ is monotone decreasing. 
Indeed, if a graph G2 is obtained by adding one edge to a graph Gi , then 

e{G2) = e(Gi) + 1, 0(^2) G {c(Gi) - l,c(Gi)}, 

so that X{G2) > X{Gi) . Hence, using induction on e(G2) — e(Gi) , 

Gi CG2^X(G2)>X(Gi). 

Furthermore 1{£(g)<l} is also monotone decreasing. (For e ^ E{G) , if e joins 
two vertices from the same component of G then £{G + e) > <?(G') obviously. If 
e joins two components. Hi and 2^2 of G, then the resulting component has an 
excess more than or equal to max{i?(ifi), 1^(7^2)} ) with equality when one of two 
components is a tree.) 

Now notice that each G on [n] is essentially a (2) -long tuple 5 of {0, l}-valued 
variables 5{i,j), ^{i,j) = 1 meaning that (i,j) € E{G) . So, a graph function f{G) 
can be unambigiously written as f{S). Importantly, a monotone decreasing (in- 
creasing) graph function is a monotone decreasing (increasing) function of the code 
S. For the random graph G{n,p), the components of S are independent random 
variavbles. According to an FKG-type inequality, see Grimmett and Stirzaker [11] 
for instance, for any two decreasing (two increasing) functions /(Y) , (7(Y) of a 
vector Y with independent components, 

E[/(Y)ff(Y)]>E[/(Y)]Eb(Y)]. 

Applying this inequality to {l/2)-^^^>l^£(^s)<L} > we obtain 



PriSn >0,£n< L) >E 



X(G(n,p))- 



E[l 



{£(G(n,p))<L}J 



Pr(5„ > 0) Pr(£:„ < L). 



n 
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Thus our next step is to determine how large £{G{n,p)) is typically, if 

(2.1.14) p=^—^ , A = o(ni/3). 



n 
For p = c/n, c < 1, it was shown in Pittel [19] that 

lim Pr(G(n,p) does not have a cycle) = (1 — c) ' exp(c/2 + c /4). 

From this result and monotonicity of £{G), it follows that, for p in (2.1.14), 

lim Pr(£:(G(n,p)) > 0) = 1. 

If A —7- — oo, then we also have 

(2.1.15) lim Pr(£:(G(n,p)) > 0) = 0, 

that is £{G{n,p)) < with high probability (whp). (The proof of (2.1.15) mimicks 
Luczak's proof [14] of an analogous property of G{n,m), with n~^'^{n/2 — m) ^ 
oo.) 

Furthermore, by Theorem 1 in [16], and monotonicity of £{G{n,p)) , it follows that 
£{G{n,p)) is bounded in probability (is Op(l), in short), if limsupA < oo. 

Finally, suppose that A — )• oo. Let C{G{n,m)) denote the total excess of the 
number of edges over the number of vertices in the complex components of G(n, m) , 
i. e. the components that are neither trees nor unicyclic. According to a limit 
theorem for C{G{n,m = {n/2){l + Xn''^/^))) from [12], C{G{n,m))/X^ -^ 2/3, 
in probability. According to Luczak [14], whp G{n,m) has exactly one complex 
component. So whp £{G{n,m)) = C{G{n,m)) , i. e. £{G{n,m))/\^ — t- 2/3 in 
probability, as well. 
Now, if 



m' = Np + 0(^^Wq). ^-=(2)' 



then 

m' = !^(1 + A'n-i/3), A' := A(l + 0{n-^/^)). 



2 

Therefore, in probability. 



£{G{n,m')) 2 



A3 3' 

as well. From a general "transfer principle" ([5], [15]) it follows then that 

£{G{n,p)) 2 
A3 3' 

in probability, too. 

This discussion justifies the following choice of L: 

0, if lim A = —00, 

(2.1.16) L = I u ^ CO however slowly, if A = 0(1), 

A^, if A ^ 00, A = o(ni/^2)_ 
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2.2 Generating functions. 

First, some basic facts about the generating functions Ceix) and Ei{x) . In- 
troduce a tree function T(x) , the exponential generating function of {A;'^~^}, the 
counts of rooted trees on [/c] , /c > 1 . It is well known that the series 

fc>i 

has convergence radius e~^ , and that 

r(x) = xe^(^\ |x|<e-^ 

in particular, T(e~^) = 1. (This last fact has a probabilistic explanation: {— rrr} 
is the distribution of a total progeny in a branching process with an immediate 
family size being Poisson (1) distributed.) T{x) is a building block for all Ci{x). 
Namely, (Moon [18], Wright [23], Bagaev [1] resp.), 

(2.2.1) C.^{x)=T{x)-\t\x), 



(2.2.2) C,{x)=\ 



1 — T(x) 2 



T^.)(6-T(.)) 

^ ' 24(1-T(x))3 

and ultimately, for all ^ > 0, 



3^+2 

(2.2.3) C,{x) = ^ ""'-^ 



i=0 



{l-T{x)f^-'^' 



Wright [23]. Needless to say, |x| < e~'^ in all the formulas. One should rightfully 
anticipate though that the behaviour of Ci{x) for x's close to e~^ is going to 
determine an asymptotic behaviour of Pr(5'n > 0,i?n < L) . And so the {d = 0)- 
term in (2.2.3) might well be the only term we would need eventually. In this 
context, it is remarkable that in a follow-up paper [24] Wright was able to show 
that 

<c fi_j^u\\M ' (q := Q,o > 0, (i£ := -Q,i > 0) (Vn>l). 

(We write X], ^jX-' <c Xlj ^j^'' when aj < bj for all j .) In the same paper he also 
demonstrated existence of a constant c > such that 

(2.2.5) c, -cQ") (£-1)!, di^c(^) £1, (£^oo). 



d=0 




£r,d 


(1- 


T{x)fr-d^ 




(6r 


- 2d)\Qd{r) 
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Later Bagaev and Dmitriev [2] showed that c = {2tt)~'^ . By now there have been 
found other proofs of this fact. See, for instance, Bender et al [3] for an asymptotic 
expansion of q due to Meerteens, and Luczak et al [16] for a rather elementary 
proof based on the behavior of the component size distribution for the critical 

G{n,m) . 

Turn to Er{x) , r > 1. It was shown in [12] that, analogously to (2.2.3), 

Er{x) 

(2.2.6) 

^r,d -25r32r-d(3^ _ ^)!(2r - d)\ ' 

where Qoir) = 1, and, for d > 0, Qd[f) is a polynomial of degree d. By Stirling's 
formula, 

(2.2.7) e, :=e,,o~(2^)-'/' ("-") r'^-'/'e-^ r ^ oo. 

Formally differentiating both sides of (2.1.8) with respect to w and equating coef- 
ficients by w^~'^ , we get a recurrence relation 

r 

(2.2.8) rEr{x) = ^^ kCk{x)Er-k{x). 

k=l 

By (2.2.3) and (2.2.6), the highest power of {1-T{x))-^ on both sides of (2.2.8) is 
3r , and equating the two coefficients we get a recurrence relation involving e^ and 

r 

(2.2.9) r£r = 2_.^^k^r-k, r>l. 

fc=l 



With these preliminaries out of the way, we turn to the formula (2.1.11) for 
Pr(5n > 0, £n 1^ L) . Notice upfront that, for L = — arising when A — )■ — oo — we 
simply have 

(2.2.10) Pr(5„>0, £:„<0) = 7V(n,p)[x"]e^(^), F(x) = ^C_i(x) + ^Co(x). 

p 2 

The next Lemma provides a counterpart of (2.1.10) and (2.2.10) for L G [l,oo). 
Lemma 2.2.1. Given L € [l,oo), 

(2.2.11) Pr(5„ >0,£„<L) = N{n,p) J] f ^ J [x"] {e^(^)F,^(x)} , 
where {F^{x)} is determined by a recurrence relation 

rAL 

(2.2.12) rF^^{x) = -J2kCk{x)F^^_,{x), r > 1, 



fc=i 
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and Fq{x) = 1 . (Here a f\h := min{a, 6} .) 
Proof of Lemma 2.2.1. Clearly 

(2.2.13) exp r- j^ w'C,{x)\ = f^ w^ F^ 

\ 1=1 J r=0 



(X), 



where F^{x) are some power series, with nonnegative coefficients, convergent for 
|x| < e~^ . This identity implies that 



Differentiating this with respect to w and replacing exp ( X]£=i ^^^^(^) ) ^n the 

left of the resulting identity with (X]^o w'^F^{x)^ , we get , after multiplying by 
w, 

(oo \ / L \ oo 

s=0 / \^=1 / r=l 

Equating the coefficients by w'' , r > 1 , of the two sides we obtain the recurrence 
(2.2.12). 

The recurrence (2.2.12) yields a very useful information about F^{x). 
Lemma 2.2.2. Let L > . For r > , 

5r f-L 

Jr,d 



(2.2.14) F^{x) = Y,- 



=0.- nx)) 



3r-d^ 



and, denoting f^ = f^Q, g^ = -f^^^ 

(2.2.15) ,, Z ..,, - r. i^^..-1 <c FH^) < 



fL L fL 

(i-r(x))3'^ (i-r(x))3'^-i -'^ '^ ^ ^-" (1-T(x))3'- 



Furthermore the leading coefficients f^ , g^ satisfy a recurrence relation 

rf\L 

2 



(2.2.16) rf^=]^Y.kckftk; /(f = l, 

fe=i 

^ r/\L ^ r/\L 

(2.2.17) rg^ = Y^kckgtk + i;Y.'''^^J'r-k^ 9o=0, 



2 A^ --»-«: 2 
fc=i fc=i 



so, in particular, f^>0 and g^ > for r > . 

Note. 1. This Lemma and its proof are similar to those for the generating 
functions Er{x) obtained in [10]. 
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Proof of Lemma 2.2.2. (a) We prove (2.2.14) by induction on r. (2.2.14) 
holds for r = as F^{x) = 1 and f^^ = f^ = 1. Further, by (2.2.12) and (2.2.3), 



FfW = ic,(,.) = 5E7I 



Cl,d 






i. e. (2.2.14) holds for r = 1 too. Assume that r > 2 and that (2.2.14) holds for 
for r' G [l,r " 1]. Then, by (2.2.12), (2.2.3) and inductive assumption. 



^ rAL 



2r 
fc=i 



rAL 3k+2 ^ 5(r-fc) „i 

</ r — k^di 



1 v^ , v^ Ck,d sr-^ 

-^2^ 2^ (i-T(x))3fe-'^ ^ (T 



,_ T(x))3fe-'^ ^^ (1 - r(x))3('^-fc)-rfi 



rAL j^L 

Ck,d J, 



'Or Z^ Z^ 



T — k^di 



2r ^ ^ (l-T(a;))3'^-('^+'^i)' 

fc=l d<3fc+2, di<5(r-fc) ^ ^ " 



Here 



< d + di < 3fc + 2 + 5(r - A;) = 5r - 2(A; - 1) < 5r, 

so (2.2.14) holds for r as well. 

(b) Plugging (2.2.14) and (2.2.3) into (2.2.12) we get 



5r fh 

Tr,d 



Sn- 



d=o ^ ^ " 



rAL , 3fc+2 5{r-k) j-L 

V"^ fC V"^ Cfc^di \~^ Jr — k,d2 



^.a 



^^ 2r ^^ (l-r(x))3fc-'ii ^^ (l-T(x))3('-'=)-'^2- 

Equating the coefficients by (1 — r(x))~'^'' (by (1 — T(x))~'^'^"'"^ resp.) on the right 
and on the left, we obtain (2.2.16) ((2.2.17) resp.). 

(c) For r = 0, (2.2.15) holds trivially. For r > 1 , inductively we have: by (2.2.4) 
(upper bound) and (2.2.12), (2.2.16), 

1 '^^'^ n f^ 

TTiL/ \ ^ ^ \ ^ 7 ^k Jr — k 

t^ [X] < ^ '" 



-Tk- 

2r ^^ 1 



r(x))3fc (i-T(x))3('-fe) 



k=l 

1 1 ''^■^ 

fl-T(x))3'-2;^ ^ fccfc/,^_fc 



fe=l 

L 



(1-T(x))3'-' 
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furthermore, by (2.2.4) (lower bound), (2.2.12) and (2.2.16)-(2.2.17), 



Pr{x)>.^^Y^k 



n 



fc=i 

rAL 



Ck 



dk 



>r 



2r ^^ (1 



(i-r(x))3fc (i-r(x))3fc-i 

Ck 



F^-k{^) 



k=l 
rf\L 



-Yk- 

2r ^ 1 



T(x))3fe 
dk 



J 1 



r — k 



9r-k 



(i-r(x))3('^-fc) (i-r(x))3('^-'=)- 






r — k 



■^^ (i-r(x))3fc-i (i-r(x))3('^-fc) 



/.^ 



(l-T(a;))3'- (l-T(a;))3'-i 



rAL 



rAL 






L 

r — k 



k=l 



k=l 



n 



9r 



(1-T(x))3'- (l-T(x))3'-i- 



To make the bound (2.2.15) work we need to have a close look at the sequence 
{fr^ar}r>o- First of all, it follows from (2.2.16) that 



f^<fr-=tT, 9r<9r:=9^- 



That is fr and —Qr are the coefficients by (1 — T(x))~3'" and (1 — T(x))~3'"+i in 
the expansion (2.2.13) for Fr{x) := F^{x). Now, using (2.2.13) for L = oo and 
(2.1.8), we see that 

2 



^r>0 ) r>0 

So, equating the coefficients by w'^ , r > , we get 



Y,Fk{x)Fr-k{x) = Er{x). 



fc=0 



Plugging (2.2.6) and (2.2.14) (with L = oo), and comparing coefficients by (1 
r(x))-3'- ((1 - r(x))-3'^+i, resp.), we obtain 



In particular. 



7 ^ fkfr-k — ^r,0] 2 ^^ fk9r-k — —^r,l- 
fe=0 fe=0 



fr < ■^£r,0) 9r < ~'^^r,l- 



RANDOM EQUATIONS 17 

Consequently, using (2.2.6) for r > 2 and d = 0, 



r — l ^ ^ r — 1 



Jr —2^r,0 — 2 2^ JkJr-k ^ o^'''" ~ 2 ^ ^^^fi-^^r-k 
fc = l fe = l 



> 



£r,0 ( 1 _ ^ V^ 

2 ~ 4 Z.^ 



'^-^M-UPS)© 



2 1' 4 ^^ V jV f^n 



>'-f{l-l/r), 
that is 

(2.2.18) ^(1 - 1/r) < ^ < ^ - ^ (1)'^'^"'^'^"''' (^ ^ °^)' 
see (2.2.7). Furthermore, using (2.2.6) for r > and d = 1 , 

(2.2.19) 5. <b ( ^ ) r'-+i/2.-'- 

And one can prove a matching lower bound for g^ . Hence, like e^ , /^ , </,. grow 
essentially as r^ , too fast for Fr{x) = F^{x) to be useful for asymptotic estimates. 
The next Lemma (last in this subsection) shows that, in a pleasing contrast, /^ , 
g^ grow much slower when r ^ L. 

Lemma 2.2.3. There exists Lq such that, for L > Lq , 



(2.2.20) /^^<, (_J, gl;<^r^—J, Vr>0. 



Proof of Lemma 2.2.3. (a) It is immediate from (2.2.18), (2.2.19) that, for 
some absolute constant A and all L > , 

fr=fr<A(^-^^\ al^ = 9r < Ar {^j-"j < r < L. 

Let us prove existence an integer L > , with a property: if for some s > L and 
all t < s , 



then 

2eJ ' ^^+1- ^ \2eJ 



/i;,<A(^V'\ ,,V<^(. + i)f^y^' 
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By (2.2.16), (2.2.21), and (2.2.5), there exists an absolute constant B > such 
that 

A function (x/L)^ attains its minimum on [0,L] at x = L/e, and it is easy to 

show that 

{e~^, X < L/e, 



Since s + 1 > L , we obtain then 



/i,<.4B(^ + ^-^;^).L- 






if we choose 



L>Li:=B'[ - + 



1 



1-e-i i_e-(3-e)/2 
Likewise, by (2.2.17), (2.2.21) and (2.2.5), 

(s+i)rf,,<^B(,+i)(f)"'L'/.x:(i)' 

so that 



if we choose 



L> L2 := {B + B'f [—^ + 



l_e-i i_e-(3-e)/2y • 

Thus, picking L = max{Li,L2} = L2 , we can accomphsh the inductive step, from 
s (> L) to s + 1, showing that, for this L, (2.2.20) holds for all t. D 

Combining (2.2.10), Lemma 2.2.1, Lemma 2.2.2, we bound Pr(5'„ > 0, <?„ < L) . 
Proposition 2.2.4. Let L € [0,oo). Then 

Si < Pr{Sn >0,£n<L)< S2. 

Here 

12;^; ^"J n 

r>0 



j:,=N{n,p)Y,I^^J[a 



(2.2.22) 



r>0 



(1-T(x))3'- (l-T(x))3'-i 
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and 

{= fr, r < L, ( = Qr, r < L, 

with fr, Qr satisfying the conditions (2.2.18)-(2.2.19). 

Note. The relations (2.2.22)-(2.2.23) indeed cover the case L = since in this 
case /o = 1 , 50 = and //" = 5^ = for r > . 

2.3 Asymptotic formula for Pr(S'„ > 0). 

The Proposition 2.2.4 makes it clear that we need to find an asymptotic formula 
for 

(2.3.1) 7V(n,p)0„,^, ^„^^:=[x"]- ^,--— , u; = 0,3,6... 

Using N{n,p) = n!g" /^(pg-s/S'jn ^^^^ Stirling's formula for n!, with some work 
we obtain 

[3 A A^ 

(2.3.2) N{n,p) = \/2^exp - ' -^^^ -^^^ 



n- +n^/^ n^/^ — 

2 2 2 



+ y + ^ + 0(n-V3(i + A4)) 

The big-Oh term here is o(l) if |A| = o(n^' "^^) , which is the condition of Theorem 
1.1. 

Turn to (l)n,w ■ Since the function in question is analytic for |x| < e~^ , 

dx. 



^^''" 2^i7rx"+i(l-T(x)) 

where F is a simple closed contour enclosing the origin and lying in the disc |x| < 
e"^ . By (2.1.10), (2.2.1)-(2.2.2), the function in (2.3.1) depends on x only through 
T(x) , which satisfies T(x) = xe^^ "> . This suggests introducing a new variable of 
integration y, such that ye~^ = x, i. e. 

k 

y = T{x)=Y,^k^-\ \x\<e-\ 
fc>i 

Picking a simple closed contour V in the y -plane such that its image under x = 
ye~^ is a simple closed contour F within the disc |x| < e~^, and using (2.2.1)- 
(2.2.2), we obtain 



•n,^ ~ i y—'e"^ exp (K{y) - f - ^ ) (1 - y)^/^- dy, 



(2.3.3) /^ 2 
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q/p ~ n, SO y~^ e^y e'^^y> would have fully accounted for asymptotic behavior of 
the integral, had it not been for the factor (1 — y)'^'^~^ . Once T' is picked, it can 
be replaced by any circular contour y = pe^^ , 9 G (— vr, vr] , p < 1. (The condition 
p < 1 is dictated by the factor (1 — y)^/^~^ .) And (2.3.3) becomes 

(/)n,w =^^(^)' 

(2.3.4) IH := r e'^^^'^) exp (-pe*74 - p^e'^'/S) (1 - pe^y/^-"" dO, 

J —IT 

h{p,0)=- fpe^^-^4— ) +npe'^ -n{lnp + i9). 

Let us choose p < 1 in such a way that, as a function of 6, \e^^^'^'\ attains its 
maximum at 9 = 0. Now |e''('''^)| = e^'^^'^^ , with 

f(p, 9) = Re h(p, 9) = - p cos 9 p^ cos 29 + np cos 6* — n In p, 

p 2p 

so that 

fg{p,9) = —p sm6'(cos6' 1. 

Then /^(p,6') > (< resp.) for 6* < (0 > resp.) if 

(2.3.5) p<^{l + np/q). 

Let us set p = e~""^ , where a = o{n^'^) , since we want p — t- 1 . Now 

-(1 + np/q) > 1 + 2"'~^^^ P < 1 - an.-^/^ + yn'^/^ 

so (2.3.5) is obviously satisfied if 

A a 

(2.3.6) ''+2-2- 

(2.3.6) is trivially met if A > 0. For A < 0, |A| = o{n^/^), (2.3.6) is met if o > |A|. 
In all cases we will assume that lim inf a > . 

Why do we want a = o{n^'^)7 Because, as a function of p, h{p,0) attains its 
minimum at np/q -^ 1 , if A < is fixed, and in this case np/q < 1 , and the 
minimum point is 1 if A > . So our p is a reasonable approximation of the saddle 
point of \h(p,9)\, dependent on A, chosen from among the feasible values, i. e. 
those strictly below 1 . Characteristically p is very close to 1 , the singular point of 
the factor (1 — y)^'^~^ , which is especially influential for large w's. Its presence 
rules out a "pain-free" application of general tools such as Watson's Lemma. 

Under (2.3.6), 

\n{p,9)\>\'/'\sm9l 
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and sign/^(p, ^) = —sign 9, so that 

f{p,0)<f{p,^)--n"^ j sin zdz 

^^•^■^^ =f{p,0)-an'/^sm\9/2) 

< f{p, 0) - aTT-^n'^/^9^ = h{p, 0) - a{Trfn'^/^9^. 

Let us break the integral I{w) in (2.3.4) into two parts, /i(w) for \9\ < 9o, and 
l2{w) for 16*1 > ^0) where 

9q = 7rn~^/^lnn. 

Since f{p,9) is decreasing with \9\ , and |1 — pe^^\ > 1 — p, it follows from (2.3.7) 
that 



(2.3.e 



/2(^)|<6(l-e-'^""'')""'e/(^'^'') 

< (ain-i/3^ -"e^^^'O) exp(-a In^ n) 

1/3/- _nr,.-l/3\ 



ai :=n^/'^(l-e-"" ' ). 
Turn to Ii{w) . This time \9\ < 9q . First, let us write 

■ a _ -1/3 

pe = e ""^ , s = a- it, t := n 



so |s| < a + vrlnn. The second (easy) exponent in the integrand of /i(u') is 
asymptotic to —3/8, or more precisely. 



(2.3.9) 



Q2{a) :-- 



1 

4 


g-.n 1 


/3 


1 

8 


e-2 


sn 


-1/3 


- 


4 


-1/3 


- 


1 

— e 
8 


-2an"^ 



=Q2(a) + 0(|t|n-i/3), 



1/3 



Determination of a usable asymptotic formula for h{p,9) is more laborious. It is 



1/3 



convenient to set q/p = ne ^"' ; thus 



^ = ^1/3 In !^ > ni/3 in(l + An'^/^) > A(l - An-i/V2), 



and 



p_A = 0(n-2/3 + n-i/3A2). 
Using the new parameters s and /i we transform the formula (2.3.4) for h{p, 9) to 



h{p, 9) = n\e 



-(At+s)n-i/3 _ lg-(^+2s)n-i/3 ^ -«„-i/3 



e "" "'"+571-1/3 



Approximating the three exponents by the 4-th degree Taylor polynomials, we 
obtain 



h{p,9) =n 



3_„-l/3A^+„-2/3/f!_^-l/f!' 

2 2 4 12 



+ Qi{p,a) 



(2.3.10) 



Qi{lJ-,a) :=n-i/3 



4! 4!2 ^ IT 

Di{t) :=n-i/3|t|(|A| + a + lnn)=^ + n-^/^{\X\ +a + Innf. 
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(Explanation: the second summand in Di{t) is the approximation error bound 
for each of the Taylor polynomials; the first summand is the common bound of 
I in + a)^ - {fi + s)^\, \{n + Is)"^ - (^ + 2a)^\ , and |s^ - a^^l , times n'^/^ .) And 
we notice immediately that both Qi{fi,a) and Di{t) are o(l) if, in addition to 
|A| = o{n^'^'^) , we require that a = o{n^'^'^) as well, a condition we assume from 
now on. Obviously 0{Di{t)) absorbs the remainder term 0{\t\n~^'^) from (2.3.9). 

Furthermore, since 



n 



-1/3 



/i = In 



np 



n-i/^A 



n 



-2/3 



A2 



n 



q I 2 

for the cubic polynomial of n~^l^ii in (2.3.10) we have 



^(y + l)+0(n-^/^(l + A^)), 



n 



3 _i/3 ^ _2/3 ^ _i^ 

2 2 4 12 



= n\- n-l^\ + nV3^ " y " ^ + 0{n-^l\X + A^)). 

Observe that the first three summands are those in the exponent of the formula 

(2.3.2) for N{n,p) times (—1).) Therefore, using (2.3.2) for N{n,p), 
(2.3.11) 

Nin,p)exp(^h{p,e)~^pe^'-^p^e^^'^ 

= [1 + 0{Dr{t))]V2^-exp (^-^ + ^ + Qip,a) + ^ + ^^ ; 
Qifi, a) -Qiip, a) + Q2{a) + 0{n-^/\l + A^)), 

and Q{p,a) = o(l) as A, a = o{v}''^'^). In particular, using (2.3.8), (2.3.10) for 
^ = 0, i. e. s = a, we see that 



(2.3.12) N{n,p)\h{w)\ <b ni/2(a^„-i/3)-«;g-ain^n ^^p ( 



\^ no? a? 



Furthermore, switching integration from 6 io t = n^'^9 , the contribution of the 
remainder term 0(Z?i(t)) to N{n,p)Ii{w) is 0{5n,w), 



n 



-1/12 



Innf'^t 



-AV6 



(ain-V3)' 



exp 



(^n 



Di{t)dt. 



(Explanation: n ^'^"^ = v}''^n ^'^n ^'^ , with n ^'^ coming from n -'^''^(a + 
vrlnn)'^/'^, an upper bound of |1 — pe'^^\^'^ , for |^| < 9q.) 

Now 






exp 



yua^ a'^ 



f+»u^ 



where, see (2.3.6), 



^ + a = -+a + 0{n-^/^ + n-^'^X^) > 0, 
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since liminfa > 0, and a > |A| if A < 0. Hence, see (2.3.10) for Di{t) , we have 
Sn,w <b ^n,w , where 



(2.3.13) A^ 



n -V12W3 . „--exp ( -^ + ^ + |- 



/3 (|A|+a + lnn)3 ^3 (|A|+a + lnn) 

/i/2 + a (^/2 + a)i/2 



• (o + lnn)^/"^ n 
The denominators ///2 + a, (^/2 + a)"*^'^ come from the integrals 



\t\'' exp 



51 



(^1 + a) t2j ^i ^ c^(^/2 + a)-('=+i)/2^ (A; > 0), 



for /c = 0, 1 . Clearly A„^^ absorbs the bound (2.3.12). 

Thus, switching from 6 to s = a — in^'^6 , it remains to evaluate sharply 



(2.3.14) 



i(2^)-i/V/6exp(-^ + ^ + Q(^,a; 



expl^ + f](l-e— -^^^)3/--d.; 



Sl 



here si = a — in^'^^Oo, S2 = a + in^'^OQ, and the integral is over the vertical line 
segment connecting si and S2 . Lastly we need to estimate an error coming from 



-1/3, 



replacing (1 — e ^" ' )^/'^ ^ with a genuinely palatable [sn i/3)3/4 w_ Using 



\sn 



-1/3, 



> ll-e 



"'""^'|>|l-e-''""^'|, {s = a-it), 



|x"-l| <n|x-l|, (u> 1, |x| < 1), 
we have: for u > 1 



1 



1 



(^l_ f^-sn-^/^Y (sn-l/3)« 



< 



1 



1 1 a — an -'-/'^ \u 



<- 



u 



n ^/^ \u 



1-e" 



,-1/3 



11 -e-"^ 



Also, for s in question. 



So 



sn" 



-1/3 



<f 



ti sn 



-1/31 



|l-e- 



-1/31 



< 



l-e" 


sn-^/^ 


Y 


sn~ 


1/3 


J 


u{a + 


\t\)n- 


1/3 



■an "-/ ■^ \u 



11 -e-"' 



1-1/3 



n ^/ -^ ti 



|l-e" 



-1/3, 



> 0.5|sn" 



-1/3, 



:i - e-^""'')3/4 _ (sn-l/3)3/4| =|^^-l/3|3/4 



1-e- 



-1/3 ■ 



3/4 



sn ^/^ / 

<,|5n-i/=^|3/4+i<,n-^/i2(^^|^|^7/4_ 
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Combining these two estimates, we have: for w € {0, 3,4,...}, 

1(1 - e-^-^^^)3/-- - (.n-V3)3/4-.| <^ (^ + l) ""y"^°v3'5^' ; 

see (2.3.8) for oi . Consequently, replacing (1 — e~*"^ ^3/4-«; jj^ (2.3.14) with 
[sn~^'^)^''^~^ incurs an additive error of order 

l/12+«,/3 . -^ ^^^f_^,[^,^\. „-l/3„5/4 



(^, + i)„-i/i2+«,/3 . „-- exp ( -:^ + ^ + y j • n-V3a 

at most; thus the error is easily 0{{w + 1)A„^^), see (2.3.13) for A„^^. 

While these bounds will suffice for A = 0(1) , the case A — > oo requires a sharper 
approximation of (1 — e~*" ^s/'i-w ^^^ ^ _ q(;\3^_ -^^e write 



(1 _ g-.n-^^^)3/4-«, ^ ^^^-1/3)3/4-^ 



exp 



(3/4 - w) In 



1 - e— ^^^^ 



sn" 



-1/3 



: (sn-i/3)3/4-«' exp[Q3(^, a) + 0{Ds{w, t))] ; 



l_g-an-/3 



(2.3.15) 

Q3(w;,a) :=(3/4-u;)ln' _ , 

an '-'■^ 

D3{w,t) :={w + l)tn-'^/^. 

Notice that Q3{a,w) -^ as wa = ©(A^n^/^^) = o{n^^^), and Ds{t,w) -> as 
tflnn = o{n^'^) . The expression (2.3.14) therefore becomes 

-i(2^)-i/2n-i/i2W3exp (-^ + ^ + Q(^,^,a)') 

■ /expf^ + ^').3/4-d, + o((u; + l)A„,^); 
(2.3.16) J V 2 3y 

Q(/x, w, a) :=Q(/U, a) + Qsiw, a); 

K^ :=n-Vi2+./3 . a- exp (-^ + ^ + ^) (a + \nnf'\ 

Finally, after this replacement we can extend the integration to (a — ioo, a + ioo) , 
since the attendant additive error is easily shown to be absorbed by {w + l)An^^ 
for all w , and by {w + l)ls.n,yj if it; = 0{\^) . 

Lemma 2.3.1. Suppose that A = ©(n^/-*^^). Let a > |A| be such that lima > 0, 
a = o{n'^'^'^) . Then, denoting fj, = n-^'"^ ln{np/q) , 



(2.3.17) N{n,p)[x 



eH{x) 



■ c2 ^3 



(1-T(a;))«' 
i(2^)-i/2e3/8+°(i)n-i/i2+-/3e-MV6 f ,3/4-^ expf^ + ^ 

a — ZOO 

+ 0{{w + l)R, 



, ds 
3 
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with Rn,w ^ ^n,w for all w , and Rn,w = ^n,w A A„^^ if wa and w\nn are both 
o{n~^'^) . Furthermore, shifting the integration line to {s = b + it: t € (—00,00)} 
does not change the value of the integral as long as b A (^/2 + b) remains positive. 

Proof of Lemma 2.3.1. We only have to explain preservation of the integral, 
and why e~^ '^ can be replaced with e~^ '^ . Given such a b, pick T > and 
introduce two horizontal line segments, Ci^2 = {s = a ± iT : a £ [a,b]} , the top 
segment and the bottom segment being respectively right and left oriented. On 

C1UC2, 

and 

- + a> - + (aA6) > 0. 

Therefore 

C1UC2 
As for e~'^ '^ ~ e~^ ' ^ , this follows from 

(2.3.18) lA^ - fi'\ <b A2(n-2/3 + n'^/^A^) = (n-i/^^)^ + n'^/^A^ ^ 0. 

n 

In the context of the critical random graph G{n, m) , the integral appearing in 
(2.3.17) was encountered and studied in [12]. Following [12], introduce 

^ a-\-ioo 

(2.3.19) A{y, ^) = ^— j s'-y exp (^^ + -J d.. 

a — icxD 

We know that this integral is well defined, and does not depend on a , if a > and 
a > — /u/2. It was shown in [12] that (1) 



3(j'+i)/3 ^ fc!r((y + 1 - 2A:)/3) ' 
(2) A(y,/i) > for y > 0, A{y,ii) > for y > 2, and (3) 

r (2^)-i/2|^|i/2-j/^ /x-^-00, 

(2.3.21) ^(y,/")~< e-^'/6 

i 2«/2r(y/2)/ii-s'/2 ^ 

We will also need two bounds 



/i — 7- 00. 



3 



^(y,^)<.e^"^lV3(2^ 
(2.3.22) Ms; 



A{y, li) <b (a + fi/2)-'^'a'-y exp 



i/2^i-?/_^/'_A^ , /f^ , ^ 
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(the second bound holding for y > 1 and a + iJi/2 > 0), and an asymptotic formula: 
if ;U — 7- oo , y — > oo , and y = 0{pt') , then 



(2.3.23) A{y, ^) ~ (27r)- ^^(y^-^ + ^ + 20"^/^^'-^ exp ( -^ + ^ + 
where ^ = ^{y, fi) is a unique positive root of 

Also, if y = O(A^), then 

(2.3.24) Aiy, /.) <, ^"V^^i-y exp (-^ + ^ + ^) , 

(See Appendix for a proof of (2.3.22) and (2.3.23)-(2.3.24).) 
With A{y,ii), we write (2.3.17) more compactly: 

(2.3.25) N{n,p)[x''' 



1-T(x))«' 
= (2^) V2e3/8+o(i)^(i/4 + ^^ /,)n-i/i2W3 + o((^ + i) A,,^). 

Let us use (2.3.25) for asymptotic evaluation of Pr(S'„ > 0, <5„ < L) given by 
(2.2.10)-(2.2.11). 

Case |A| = 0{1) . According to (2.1.16), we can pick L — )• oo as slowly as we 
wish. We pick L = ln ' n. 

As a first step, let us estimate the overall contributions, Rn and Rn , of the 
remainders 0{{w + l)Rn,w) to the bounds Si and S2 in Proposition 2.2.3. In this 
case we choose a = {Ly^^ for each w , and Rn,w = A„^^ . Consider Rn first. By 
(2.2.19) and (2.3.13), and dropping {3r + l){np/2qy = (3r + l)(l/2 + o(l))'' factor, 

(2.3.26) __ f A3 , ^a2 , a^ 

r=0 



■«Pi-^ + ^ + y)-E/^r''- 



00 



Now ai ~ a, so by (2.2.23) and (2.2.18), 

E/.V''s.E(5i?)' + E(|^)'sE(!i <- 

r=0 r<L ^ ^' r>L ^ ^' r>0 ^ ' 

So (2.3.26) becomes 

4') <b n-1/12 . ^-l/3gln3/^n ^ ^-5/12+0(1) _ 
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Further, by (2.2.22) and 5^ = 0, 



.,_.,.^,,„,,)g(|)'[.»,_ 



gL^H{x) 



Therefore 



V- r(x))3'-i' 



.15/^ n) 



So, using the bounds (2.2.19) and (2.2.23) for ff^, we conclude that \R^n\ < 21^^ . 
Thus, for L = ln ' n, 

(2.3.27) St + 0(n-V3+o(i)) < ^;(^^^^> 3";3^"_^^^f ^ < E^ + 0(n-V3+°(^)). 



where 



(2.3.28) 



^2 = Ef?^)'/''^(V4 + 3r,/i), 

S* =S* - n- V3 j; (^) ' 5.^^(-3/4 + 3r, ^). 



r>0 

Let us have a close look at S* and Sj . Write 



^2 = E (f ) ' ^^( V4 + 3r, /.) + j; (^) ' /.M(l/4 + 3r, /.) 

— Zj21 -t ^22- 

By //- < /r , (2.2.18), (2.3.22), and Stirling's formula for T{r) = (r - 1)!, 

2^2 <.e^l'^l^/^E Q + 0(|A|--^/^))' 0)^'^-^/^e- 0)'r-i(r) 

Further, since uniformly for r < L, 

(1 + Xn-^/^Y = exp(0(L|A|n-i/3)) = 1 + 0(n-V3 jj^i/4 ^^^ 

we have 

S;i = (l + o(l))J]^^(l/4 + 3r,^). 

r<L 

And, analogously to S22 ' 

5]^A(l/4 + 3r,^)<6(^)^ 



Therefore 

^2 - E i^^(l/4 + 3r, ^) ^ E i^^( 1/4 + 3r, /x). 



^^(l/4 + 3r,^)^5].^'- 

r<L r<L 



Also, by the definition of Si in (2.3.25), it follows that \T.{ - S^| is 0{n-^/^). 
Hence 
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Proposition 2.3.2. For \\\ = 0(1), 

(2^)1/2,3/8^-1/12 -<^) :=^-^(l/4 + 3r,^), 

and /i (= \ + 0{n~^'^)) can be replaced with \, as c{x) is positive and continuous 
for all X. 

Case A ^ oo, A = o{n^^^^) . According to (2.1.16), we select L = aX^ , a > 
2/3. This time we use a refined version of (2.3.24), with the exponential factor 
sneaking behind the sum operation for r < aA^ , which allows us to choose a (< 2 A) 
dependent on r for r < aA'^ . Also, for those r and a, ra = 0{X'^) = o[n^''^) and 
rlnn = O(A'^lnn) = o{n^'^); so Rn,3r = ^n,3r in this range. For r > aX^ we 
select a = A, and here Rn,3r = ^n,3r ■ (So, a = o{n^^^'^) throughout.) By (2.2.20), 

^"p(-Y + ^ + y) E(- + i)(V2+o(i))v.^ 

/2A3\ ^ /3«A3(l + o(l))y 



^ .3 /2A3 3a(l + o(l)) 

<feA'^ exp h A'^a In — ^^ ^^^ 

\ 3 4e 

. ,3 ^3 1 3a(l+o(l)) 
<bA exp I A a m 

and, pushing a down to 2/3, we can make the coefficient of A'^ in the exponent 
arbitrarily close to 

2 , 1 

- . In - = -0.46 .... 

3 2 

According to (2.2.18) and (2.3.13), it remains to bound 

(we have replaced Oi = n^'^{l — e~°'"' ) with a, since for r < aA'^ , 

and A'^/6 with /i'^/6, see (2.3.18). So we need to find minQ<2A ^{r,a) , 

(2.3.29) c|>(,,a):=-^-rlnf^V — + -• 
^ ^ ^ ' ^ 6 V4eaV 2 3 

^{r,a) attains its absolute minimum at ^(r), a unique positive root of 

(2.3.30) ^^ + ^2^ < 
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i. e. ^(r) < 2A if a is sufficiently close to 2/3 from above. Further (j)(r) := 
<l>(r, ^(r)) attains its maximum at f, a root of 

0'(y) = ln^-31n^(y)=O, 



1. e. 



4 .,3 . . ^..x /^ 



(2.3.31) ^=^^' «:=e(r) = f (<2A). 

Consequently 

(2^32) ,„).exp(-i5^)~e.p(-15^) 

It is easy to show that 

27 27 






(2.3.33) 0"(r-)- ^3 - ^3. 

and, with some work, that (j)"{r) < always. A standard application of Laplace 
method yields 

Therefore, consulting (2.3.16) for An,w and (2.3.13) for A^^^ , we bound i?„ , the 
total contribution of the remainders {w + l)Rn,w to the sum E2 in (2.2.20): 

R(^) = Y^ir + l)fi„,3. <. n-Vi2 exp (-^) {n-'/'\'-'' + n'^s in^/^ ^) 

(2.3.34) <b (A-i/4 + n-i/3+°(i))n-i/i2 exp (-^) • 

As for Rn , the total contribution of the remainders {w + l)Rn^w to Si in (2.2.20), 
it is 0{Rn ), just like the A = 0(1) case. So we arrive at the counterpart of 
(2.3.27)-(2.3.28), with 



(A-V4 + n-V3+o(l))^^p(^_10A3^| 



taking place of n i/3+o(i) . Further, again we split 112 = S21 + S22 • To bound II22 
we use the second bound for ^(1/4 + 3r, /i) with a = X (2.3.21), and the bound 
(2.2.20) for /^. Just like R^n\ we obtain 

2^2 <b A^/2 exp (-— + — + — 
(2.3.35) • J]] (l/2 + o(l))7,^<bA'/^e-°-^^^' 
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To evaluate sharply S21 , we use (2.3.23) to approximate A{y,fi) for efx < y, 
y = 0{\^), and (2.3.24) to bound A{y,^) for y < efi, e > sufficiently small. 
Invoking (2.2.18) as well, we have 



eM<T^aA3V ^1 J '^^e^<7^„A3 (r((3r + l)^-^ + ^ + 2^) 



1/2 



here 0(r) := mma^{r,a) = ^{r,^), see (2.3.28)-(2.3.29) for ^{r,a) and ^ = ^(r). 
We know that (/)(r) attains its pronounced maximum at f = (4/27)A^, i. e. well 
within [efi,aX^]. Using (2.3.31)-(2.3.33), by Laplace method, 

V- fnp\\ ^, , ,1 i^l^ ( 27r ^'/' 



£/j,<r<Q:A 



1 1/4 / 10A3 

A ' exp 



4(27r)i/233/4 ^ ^ 81 

Applying (2.3.24), it is not difficult to show that 

So 

^-= E(i)^-^(V4 + 3,...)~jp^ii^AV^exp^ >°^' 

r<QA3 ^ ^ / ^ ^ 

hence (see (2.3.35)) 



4(27r)i/233/4- ^\^ 81 

as well. And, analogously to the A = 0(1) case, for S* defined in (2.3.28), 

(2.3.37) |St - S;| « A^/^exp (-^) • 

Proposition 2.3.3. For A ^ oo, A = o(n^/-'^^), 

p3/8 / inx3\ 

Pr(5„ > 0, £:„ < L) ~ ^ AV^xp ^— ^j . 



Proof of Proposition 2.3.3 The probability is asymptotic to the expression 
in (2.3.36) times (27r)i/2g3/8^-i/i2 _ q 

Lastly, 



RANDOM EQUATIONS 31 

Case A — t- — oo, |A| = o(n^'^^). According to (2.1.16), we can pick L = 0. By 
Proposition 2.2.4 and (2.3.17) for w = 0, and a > |A| , a = o{n^''^'^) , we have 

Pr(5„ > 0,£:„ < 0) =N{n,p) [x"] e^(") 

= i2n)'/'n-'/''e'/' Ail/4, f,) + OiA^,o). 



Notice that 

' ua^ a^\ A^ 



a=|A| 6 



since X^ — fi^ = o(l). Setting a = A in (2.3.13), we obtain 

A.,0 « n-l/12 . ^-1/3 . ^-l/3(|^|3.75 ^ 1^3.75 ^y 

And, by (2.3.21), 

A{l/4,fi) ~ (27r)-i/Vl'/'"'/' - (27r)-i/2|A|i/^ 
Proposition 2.3.4. Suppose A — > — oo, |A| = o(n-^/^^). T/ien 

Pr(5„ >0,£:„<0)~e3/8|A|i/l 

Since in each of the three cases our L is such that 

Urn Pr(£:„ < L) = 1, 

Propositions 2.3.2 -2.3.4 combined with the relations (2.1.12) and (2.1.13), prove 
the part of Theorem 1.1 about G{n,p), p= 1/2. 

3. Solvability probability: G{n,m) and p = 1/2. 

Our task is to show that the result for the near-critical G{n,p), p = (1 + 
Xn~^'^)/n, A = 0(77."^'^^), implies the analogous claim for G{n,m), m = (n/2)(l-|- 
An-i/3). Denoting A^ = Q) , 

(3.1) p=!^+o{m'/^N-^) = ^-±^^^^, X' = X + 0{n-'/^)). 

N n 

Obviously A' = o{n^'^'^) , so 

Pi{S{G{n,p) >0) ^0. 

Since an event {S{G) > 0} is monotone (increasing) with G, a general "p-to-?Ti" 
result, Bollobas [5], Luczak [15], implies that 

Fr{S{G{n,m) > 0) -> 0, 
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too. However we want to prove a sharp formula 
(3.2) Fi{S{G{n, m) > 0) ~ c{X)n-^/^^, 

so that the probabihties in question can be as small as 

exp ('-l^(nVi2-o(i))3A ^ g^p j'_l^„i/4-o(i) 

It turns out that in our case the argument in [5], [15] can be sharpened to yield 

(3.2). 

To start, recall the classic entropy bound 

Pr(Bin(7V,p) > k) < exp[NH{k/N)], k > Np, 
Fi(Bm{N,p) <k)< exp[NH{k/N)], k < Np, 

where 

H{x) : X ln{p/x) + (1 — x) ln(g/(l — x)). 

Approximating H{x) by its second degree Taylor polynomial plus a remainder 
term, we obtain: uniformly for p < 1/2, and uj < a{NpY'^ , a > being fixed, 

,2 , 



Pr (Bin(iV,p) >Np + uj^/Npq] <b e""^ /^ 

(3.3) ; ; 

Pr (Bin(iV,p) < Np - uj^/Npq] <be-"'/^ 



(The bounded factor implicit in <b notation depends on a.) Given m and oj < 
m^/^ , introduce pi < P2'- 



(3.4) 
Then 



Npi +uj-\jNpi =m =^pi = (4iV)"^(-v/4m. + tj2 ~ujf, 
Np2 - uj\/Np2 =m ^ p2 = {4:N)~^ [\/4rn+^ + ujf . 



Np2 Npi m 



(3.5) i^" W 



6 



UJ^ 



(^/TTujyAm - uj/{2^/m) 



>a:=(V2-l)2, 



as u)/2^/m < O.bm-^/^ < 1 . 

Now, using e(G) to denote the number of edges in a graph G, e{G{n,p)) 
Bin(iV,p). So, by (3.3)-(3.5), 



Pr(e(G(n,pi)) > m) < Pr(e(G(n,pi)) > Npi+u^/N^i) <b e""^ /^, 

(3.6) 

Pr(e(G(n,p2)) < m) < Pr(e(G(n,p2)) < Np2 - uj^/N^) <i, e''^ /^. 

Since 

N 

PT{S{G{n,p)) >0) = Y^ Pr(e(G(n,p)) = ^) PT{S{G{n, fi)) > 0), 

11=0 
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and Pr(5'(G(n, /i)) > 0) decreases with /i, we have 

Pv{S{G{n,pi)) > 0) > Pr(e(G(n,pi)) < m) Pr(5(G(n,m)) > 0) 
>(1 - 0(e-^'/2)) Pr(5(G(n,m)) > 0), 

and 

PT{S{G{n,p2)) = 0) > Pr(e(G(n,p2)) > m) Pr{S{G{n,m)) = 0)) 
=(1 - 0(e-^'/2)) Pr(5(G(n,m)) = 0). 

Therefore 

Pr(S(G(n,j92)) >0)-O(e-"'/') 
1 - 0(e-^'/2) 

Now, by (3.4), 

Pi,2=^(l + 0(^rn /)) = ^^ , 

Ai,2 =A + O{ojm-^/^ + n-2/3), 

so, as |A| = o(n^/^2), 

\l^ =X^ + [X'iivm-'/' + n-2/3)] + Q [{um-'/'f + n'^] = A^ + o(a;n-i/3)_ 
That is, Af 2 — A^ — > 0. Hence, 
(3.8) Pr(5(G(n,pi,2)) > 0) - c{X)n-^/^\ 



Also u;^ ^ |A|'^ if u; = n^'*, which is compatible with the restriction to < n 
For this choice of u , the relations (3.7)-(3.8) imply: for A = o{n^'^^), 

Pi{S{G{n,m)) > 0) ~ c{X)n-^/^^. 

This completes the proof of Theorem 1.1 for p = 1/2. D 

4. Solvability (2-colorability) probability: G{n,p), G{n,m) and p=l. 
Consider the G{n,p) case. We know that the system 

Xi + Xj = l(mod 2), {i,j)eE{G) 

is solvable iff the graph G has no odd cycles. So a counterpart of (2.1.11) is 

" L 
(4.1) Pi{Sn >0,£n<L)= N{n,p) [x"] exp 



1/6 



E (f) c?(.) 

£=-1 ^^^ 
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{Sn = S{G{n,p)) , £n = £{G{n,p))), where C|(x) is the exponential generating 
function of graphs G without odd cycles, with an excess £{G) = i. And again the 
events {Sn > 0} and {<?„ < L} are positively correlated, i. e. 

Pr(5„ > 0, £:„ < L) < Pr(5„ > 0) < ^'^^" f °'/" " ^^ 

Thus the generating functions C|(x) take a center stage. Obviously 

Furthermore, while 

for Cq{x) we have 

(4.2) ^o(^) = iO"T^^"^'^"^ 

Indeed, we enumerate the connected unicyclic graphs with an even cycle, i. e. 
forests of an even number of rooted trees, whose roots form an undirected cycle. 
So 

(j-l)!T^(x) 



C^{x)= Y. 



2 j! ' 

even j/'>4 

which simplifies to (4.2). Comparing Cq(c) and Cq{x) we see that, for |x| < e~^, 
X — 7- 1 , i. e. for x dominant asymptotically, 

(4.3) C^{x) = lco(x) + ^ ~ J ln2 + 0(|T(x) - 1|); 

in particular, Cq{x) ~ (l/2)(7o(x) . We want to show that this pattern persists for 
^ > 0, namely 

(4.4) C!{x)^^Ce{x), {\x\<e-\x^e-'). 

Comparing (2.1.11) and (4.1), and recalling the different roles played by Cq{x) and 
{Ce{x)}iyQ in the analysis of the p = 1/2 makes it transparent, hopefully, that for 
p = 1 we should have 

Pr(G(n,p) is 2-colorable) = Pr(5„ > 0) ~ 2-^/^e^/^c{X)n-^/^^. 
Let us prove (4.4). First 
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Proposition 4.1. Given n and m < N := (2)? ^^^ C{n,m) denote the total 
number of connected graphs on [n] with m edges, and let C'^{n,m) denote the total 
number of connected graphs without odd cycles. Then 



(4.5) 

Consequently 

(4.6) 



1 



C^in.m) < — Cin.m). 



Cl{x)<c,^,C,{x), l>-l. 



Proof of Proposition 4.1. We begin with a simple claim. 

Lemma 4.2. Let T be a tree on the vertex set [n] . Let X{T) denote the total 
number of paths in T of an even edge-length 2 at least. Then X{T) > X(P„), 
where Pn is a path on [n] , and 



(4.7) 



X(P„ 



n{n — 2) 



Proof of Lemma 4.2. Pick a vertex u G [n], and introduce Vo(T) and Vi{T) 
the set of vertices reachable from v by paths of even length 2 at least, and odd 
length respectively; in particular v £ Vq. Now every two vertices from Vi{T) , 
{i = 0,1), are connected by an even path, while there is no even path connecting 
vo G Vo{T) and vi G Vi{T). Hence 



X{T) 



Vo{T)\ 
2 



Vi{T)\ 
2 



It follows that X(T) attains its minimum when |Vo(T)| = \_n/2\ and |Vi(T)| 
\n/2\ , or the other way around, i. e. when T = Pn, and the minimum value is 



XiPr, 



[n/2\ 
2 



[n/21 
2 



n{n — 2) 



n 



Armed with this Lemma, we will derive a recurrence inequality for C^{n,m). 
First we recall a recurrence equality for C{n,m) , [23], [3]: for n > 3, n—l<m< 

N, 



(4.8) mC{n,m) = {N -m+l)C{n,m-l) 

1 v^ 



+ 



2 ^^ 

m-^ -\-m'2^Tn — 1 



n 



Ui 



nin2C{ni,mi)C{n2,m2) 



Explanation. The left hand side of (4.8) is the total number of the connected (n, m) 
graphs with a marked edge. Each one of these graphs with a marked edge can be 
obtained in one of two, mutually exclusive ways. First way is inserting a marked 
edge into a connected graph on [n] with m — 1 edges, which accounts for the first 
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term on the right hand side of (4.6); indeed N — m + 1 is the total number of 
unordered pairs of vertices not connected by an edge in a given connected graph 
with m — \ edges. Second way is to start with a connected (ni,mi) graph and a 
connected (n2,?Ti2) graph, having m — 1 edges in total, and to add a marked edge 
that joins two connected graphs; nin2 is the total number of ways to select two 
"contact" points, representing each of two graphs. 

Let us see if there is a similar recursive formula for C'^{n, m) . Clearly, if a marked 
edge joins two connected graphs, none of these two graphs may have an odd cycle. 
So we definitely have the "C"^(-,-)" counterpart of the second term on the right 
hand side of (4.8). As for a potential counterpart of the first term, a difficulty 
is that an additional ?Ti,-th edge is not allowed to form an odd cycle with any of 
the m—1 edges already present. And so the total number of admissible options 
depends on the structure of a (n, m—1) graph G in question. (For such a graph to 
be connected, it is necessary that m'>n.) However we can bound the number of 
options. G is spanned by a tree T on [n] , and none of the ?ti— 1 — (n — 1) = m — n 
edges of G\T completes an odd cycle by joining the ends of an even path in T. By 
Lemma 4.2, the total number of those even paths is \n{n — 2)/4] , at least. Hence 
the total number of options for the m-edge \s N — (m — 1) — \n{n — 2)/4] , at most. 
And it is straightforward that, for n>3 and by m,> n, 



N -im-l] 



nin 



<i(Ar-(m-l)). 



So C"^(-, •) satisfies a recursive inequality: for n>3, n — 1<?tt,<A'^, 
(4.9) mG%n,m)<-{N-m + l)G^{n,m-l) 



+ 2 ^ ( )^i^2C"'(ni, 771,1)^(712, m2). 



"•1 +''T'2 — "-' 



{G''{i',fi) := if z^ = 0, or // ^ [z/ - 1, (2)]-) We will use (4.9) and induction 
to prove the bound (4.5). To this end, we define a lexicographical order, -<, on 
{(77,777) : 77 > 1, 77 — 1 < 777 < (2)} as follows: denoting £ = m — n, 

(77i,777i) ~< (772,7772) <=^ ii < £2, or ii = £2 and 77i < 772. 

The order -< is total, and (1,0) is the minimal element. The inductive basis holds, 
since C^(1,0) = C(1,0) = 1 , and C^(2,l) = C(2, 1) = 1. Suppose that, for some 
77 > 2 and 777 G [77 — 1, (2)] , 

C"'(^' m) < ^^_^+i C{u, fi), V {u, n) -< (77, 777). 
Since (77, 777 — 1) ^ (77, 777) , the inductive assumption implies that 

^(^N-m + l)G%n,m-l)<^(N-m + l) ^^_l_^^^ Gin,m-l) 
(4.10) =_±_^N-m + l)G{n,m). 
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Further, for the double sum in (4.9), 

"T-i — JT-i + 1 > 0, m2 — n2 + 1 > 0, 
and 
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SO that 



(mi — ni + 1) + (?Ti2 — n2 + l) = m — 1 — n + 2 = ?7i — n + 1, 



rui — Ui + 1 < m — n + 1 =^ rui — rii < m — n, i = 1,2. 



So, for m, 77-2 > 0, we have {ni,mi) -< {n,m) and therefore, by the inductive 
assumption, 

(4.11) llC^{n,,m,) < n^;^II^C'(n„m,) = ^^^U^i^^^^^)- 

i = l 2=1 i=l 

Combining (4.9)-(4.11), and the recurrence equation (4.8) for C{-, •) , we obtain 

1 



mC%n, m) < „ , ^ (N - m + l)C{n, m - 1) 



Om — n-\-l 

1 1 
+ 



^ (:)n-^( 



2m-n+l 2 ^^ ym, 

TTij^ -\-m2^'m — l 

1 



ni,mi 



2»Ti — n+l ^ ' ' 

Thus the bound (4.3) holds for {n,m) too. The proof of Proposition 4.1 is com- 
plete. D 

By Proposition 4.1 and and the formula (4.1), we have 



(4.12) FriSn >0,£n< L) 

< N{n,p) [x"'] exp 



^Ci(x) + Conx) + iX;(|^)'^K= 



Since Cq{x) is asymptotic to (l/2)Co(x) + ln(2 ^/^e^'^) as x — )■ e ^ , only a trivial 
change in the proof of Theorem 1.1 (i) is needed to show that 



(4.13) 



Pr(5„ > 0) < 2-i/^ei/«c(A)n-Vi^ (|A| = oin'/'^). 



We omit the details. Furthermore, since for A — t- — oo we use L = 0, the sums 
^g^i in (4.1), (4.12) disappear, and we obtain an asymptotic equality 

FiiSn > 0) ~ 2-i/4e^/*c(A)n-i/^2^ (|A| = o{n'/^\ A ^ -oo). 

To complete the proof of (ii), (case A = 0(1)), we need to prove (4.4) for each 
fixed £ > 0. Recall Wright's formula 



(4.14) Q(x) = (l-r(x)) 



-31 



21 



^c,Ai-Tix)Y 



.d=0 



, (^>o). 



Let us find a similar formula for C^{x) , £ > 0. 
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Proposition 4.3. For £ > 0, 



2r^\\-3e 



(4.15) Qix) = (1 - T\x)) 



^4Al~T{x)Y 



d=0 



where 

(4.16) c,% = 22^-1 Q,o. 
Consequently, for \x\ < e~^ and x — t- e~^ , 

C!{x) ~ ^Ce{x). 

Proof of Proposition 4.3. We use the ideas of Wright's original proof of 
(4.14), and the improvements suggested by Stepanov [22], (cf. [12], Section 9). 

Given a connected graph G on [n] , with an excess i = e{G) — v{G) > 0, 
we apply a "pruning" algorithm which successively deletes vertices of degree 1 . 
Obviously the excess is preserved, and so for a terminal graph (core) G we have 
e{G) — v{G) = i. G inherits all the cycles of G , and thus G has only even cycles 
iff G does. A minimum degree of (5 is 2 at least, and — since i{G) = i > — a 
maximum degree is 3 at least. Next we apply a "cancellation" algorithm to G: at 
each step, we delete a vertex of degree 2 , splicing together the two edges it formerly 
touched. The excess is preserved again. Once all the vertices of degree 2 are gone, 
we get a connected multigraph (kernel) G, with possible loops and parallel edges, 
and a minimum vertex degree 3 at least. Thus 

2e(G) > 3t;(G'), e{G) - v{G) = i, 

and so 

(4.17) v{G) < 2i, e{G) < U. 

Notice that the largest numbers of vertices and the edges in the kernel are 21 and 
3£ respectively, and the corresponding kernel is a 3 -regular multigraph. (In [12] 
graphs G with such kernels were called clean. It is these clean graphs that are 
most populous asymptotically among all connected graphs on [n] with excess £.) 
Now that we have a reduced number w(G) of vertices, we relabel them using indices 
from [11(0)] and preserving the order of their old indices from [n] . Under this rule, 
it follows from (11) that the number of kernels G for the collection of all connected 
graphs G on [n] with excess £ is a function of i only! 

A key element of Wright's argument was the following identity. Let M be a 
connected multigraph on a vertex set \u\ , with [li indistinguishable loops at vertex 
i, and yuij indistinguishable parallel edges joining i and j, (i,j G \v\^ i 7^ j). 
Let hn^M denote the total number of the connected simple graphs G on [n] , with 
minimum degree 2 at least and maximum degree 3 at least (core-type graphs, in 
short), such that G = M . Letting 

(4.18) 7^^(,) = ^Vm^.^ 
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we have 



(4.19) Um^z) = — •/Cm(^), ^:= J]]Aii + Xl^*J' 



where 



(4.20) 



I/! (1 - zY 



icm{z)= n ^"'^ n ^^'^"'(M^i-(/"^i-i)^) 






(Observe that /Cm(1) = 1-) Once (4.19) is estabhshed, it is easy to determine 
Hm{x) , the exponential generating function of all connected graphs G whose kernel 
is the multigraph M . Indeed to go from a core G back to G on [n] we need to 
choose an ordered sequence of v{G) of rooted trees, of total size n, and plant them 
at the vertices of G , moving increasingly from vertex 1 to vertex v{G) . Since the 
generating function of such sequences is T{xY^'^' , we see that 

(4.21) Hm{x) = V ^ T(x)" = UMinx)). 



Finally 



(4.22) 



Y,cU-^-T{x)y 



.d=0 



2e 
C,ix)= Y. Hm{x) = {1-T{x))-^' 

M:e{M)-v(M)=e 
^ ^ II meets (4.18) 

Our first step is to obtain a counterpart of (4.18)-(4.20) for 

n 

where h^ j^ is the total number of the connected core-type graphs on [n] with 
only even cycles, that cancel to a given multigraph M . To this end, consider an 
auxilliary problem. Let 

(4.23) ^i, = ^i<l + ^i°, fn, = ij.'ij + fi"^, {i<i^j<u). 

Let /in,(^^,/x°) denote the total number of the core-type graphs G on [n], which 
cancel to M , such that: (1) for each i, G has an even (odd resp.) number of 
2-degree vertices put on each of fif (/U° resp.) loops at vertex i of M; (2) for each 
(i,j), G has an even (odd resp.) number of 2-degree vertices put on each of //?■ 
(/U°- resp.) parallel edges joining the vertices i and j in M . Let us determine 
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A core-type graph G on [n] cancelling to M and meeting the parity conditions 
(l)-(2) can be viewed as a partition of [n] into: 

(a) a subset of cardinality v , whose elements are the assigned to the v vertices of 
M in a unique (order-preserving) fashion; 

(b) Vi e [v], a collection of fii ordered subsets, each having 2 elements at least 
(as G is simple), such that exactly /^| {fi° resp.) subsets have an even (odd resp.) 
number of elements; 

(c) V 1 < z 7^ j < 1/ , a collection of fiij ordered subsets, with at most one empty 
subset (as G is simple), such that exactly /if- (/i°- resp.) subsets have an even 
(odd resp.) number of elements. 

So 'H(^e^^o)(z) is the product of generating functions 'Hf(z) corresponding to 
1 + z/ + (2) combinatorial structures described in (a), (b), (c). The first is easy: 



(4.24) 

Next, for i € [z^] 



n{z) = -. 



nAz) 



2^-'^xl\^^1\ 



E' 



^",(m-.m°)^ ! 



here an,(/i'=,^°) is the total number of compositions of n with iXi parts, each 2 at 
least, such that the first /if parts (the last /U° parts resp.) are even (odd resp.). 
The factor 1/2^^ is needed as we do not distinguish between two opposite orderings 
of vertices sprinkled on each of fii loops of M at i. Consequently 



n,{z) 



(4.25) 



2A'«/if!/if! 



1 




2A''/if !/i°! (1 - z2)M« 
Similarly, for each 1 < i < j < i^ , 

1 



(4.26) n^Az 



,,6 \,,0 I 



\~Z^ 



^ 



2fc + l 



«J 



1-Z2 



M^,-l' 



1-Z2 



A'r, 



Taking the product of the generating functions in (4.24)-(4.26) we obtain 



(4.27) 
where 



'^(m%m°)(^) 



V\ (1-Z2)M 



•^(^%^°)(^), 



(4.28) /C(^.,^o)(z) 



n V - 



^A^2a*.+m° . 



n 



Ml 



\<i<j<v ^^»J 



z'^-^+^^-2[/i^-(/.^-l)z2]. 
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As a partial check, summing over (/i'^,/j°), we obtain Wright's formula (4.18)- 
(4.20). 

Now, for a core-type graph G on [n] without odd cycles, that cancels to M, 
G's parity parameters /2^, /x° must satisfy certain conditions. First of all, for each 
i G [z^] , //? = 0, since otherwise G would have an odd cycle, with a single branching 
vertex. Likewise, for ^ij > 0, the numbers of 2-degree vertices of G on fiij parallel 
edges of M must all be of the same parity, hence /uf- = ^Xij or ^°- = jiij . Subject 
to this condition, how many choices for (/j'^,/j°) do we have? For each {i,j) such 
that Hij > 0, define 

_ r 1, if ^'j- = Aiij, 

If C is a cycle in M , then the parity of a cycle in G that cancels to G is the parity 
of 6(C) := ^/j ,)gc'6ij- Hence 6(C) must be even for all cycles C, and we need 
to check this condition only for simple cycles that do not use parallel edges. Let 
T = T{M) be a tree on [i^] that spans M . Pick ^f- for all z/ — 1 pairs {i,j) such 
that {i,j) G E{T) , i. e. one of fj,ij parallel edges is in E(T). Let fj,ij > and 
6 = (^) j) ^ E!{T) . Then e completes a cycle C with a path in T that connects i 
and j . The condition "6(C) is even" determines ^f- uniquely. Hence a choice of 
z/ — 1 values of ^f ■ determines uniquely the remaining fj,'^, . Arguing as in the proof 
of Lemma 2.1.1, we see that the condition ^^h{G) is even" will hold for all other 
cycles C. Thus we have have 2^~^ choices for (fi^,fi°). 

For each of those choices, (4.28) becomes 

(4.29) /C(^.,^o)(z) = n ^"^' • n ^'^''''~'\^^^J - i^^^J - 1)^'] • n ^^"■- 

For each of these 2^~^ polynomials, 

(4.30) /C(^e,^o)(l) = l, 

and 

deg/C(^e^^o)(z) < 3^^i + 2 ^ ^Xij. 

Using the constraints 

we easily obtain then that 

deg/C(^e ,^o)(z) <Aix- ?>u. 

Now, the generating function T-L\j{z) of the core-type graphs G without odd 
cycles that cancel to M is the sum of 'H(^e^^o)(2:) over all 2*^"^ sets of feasible 
pairs (fi'^,fj,"). Using (4.27), (4.29) and (4.30) we arrive at the following formula. 
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Lemma 4.4. For each kernel M , with i := fi — i^ > , 



(4.31^ 



T-1-m{z) 



-Pm{z), 



where Pm{z) is a polynomial of degree Afi — Su = ji + 3i at most, and Pm(1) = 1 • 

This Lemma directly implies 
Corollary 4.5. 

(4.32) 



Q(x) = Yl ^m(^), 

M:e{M)-v(M)=e 



where 
(4.33) 



w-i 



Hl,{x) 



u\ {l-T'^{x)Y 



Pm{T{x)), 



Using e < fi{M) < 31, u{M) = ^(M) - £, we deduce from (4.32)-(4.33) that 

"8^-1 



(4.34) 

where 

(4.35) 



C!{x) = {l-T'{x)) 



-3i 



Y4Ai-T{x)Y 



d=0 



>2^-l 



C^,0 



i2i)\ 



Y^ K(/i). 



fj, meets (4.19) 



So, by the second line in (4.22), c^q = 2^^ '^ci^. The proof of Proposition 4.3 is 
complete. D 

Comparing (4.14) and (4.34)-(4.35), and using T{e~'^) = 1, we obtain: for i > 0, 



(4.36) Q{x) = ^^C,{x) + 0{\l-T 

And we recall, (4.3), that 



-3^+1^ 



Ixl < e ^, X — )• e "'") 



1 



(4.37) 

Now, by (4.1), for a fixed L > 0, 



Coix) = -Co{x) + ln(2-V4ei/8) + o{\T{x) - 1|). 



Pr(5„ >0,£n<L) = N{n,p) i x"""! exp V (?-\ C\ 



\x\ 



dx, 



where F is within the disc |x| < e "^ . As in Section 2.3, we switch to y by x = ye 



and choose in the y-plane the circular contour T' y = e 
fixed this time. Observe that, for each 1 < i < L, 



-i/3+ie 



o > being 



|1 - T{ye-y)\-''+' <, n-'\l - y\-''+' <, n-^'\ 
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and, likewise in (4.3) the remainder term 0(|T(ye~^) — 1|) is 0{n~^'^). And of 
course C^i{ye~y) = C_i(ye~^). On the basis of (4.36)-(4.37), it can be shown 
then that 

Pr(5„ >0,£n< L) 



2-VV/8iV(n,p)jf (ye-s')— lexp if] (^^) C,{ye-y) 



t=-\ 



d{ye~ 



where now V can be replaced by a circular contour of an arbitrarily small radius. 
Going back to the x-plane, we recognize (see (2.1.11)) the value of the resulting 
integral as 

2-1/ V/« Fl{Sn>0,£n<L)\-^i/2, 

for A = 0(1) needless to say. By (2.1.13) the latter probability is at least 



PriSn > 0)1^^1/2 • Pr(£:n < L). 
Letting n — )• oo, and using the part (i) for Pr(S'„ > 0)\p^i/2 , we get 

Since Pr(S'„ > 0) > Pr(S'„ > 0, <?„ < L) , and <?„ = Op(l) , letting L f oo enables 
us to conclude that 

Fi{Sn > 0) > 2-i/4ei/«c(A)n-i/^2_ 
Together with (4.14) this proves that 

PiiSn > 0) ~ 2-i/4ei/«c(A)n-i/i2_ 

The proof of Theorem 1.1 (ii) is now complete. D 
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Appendix. 

Proof of (2.3.22). (i) For the second bound, we use (2.3.19) and, setting s = 
a-\-it, 



exp (^ + y ) - "'"'' ^""p irr + y ) ^^p[-*'(« + /^/2)] • 



So 



A{y,^i)<{2^)-\[^ 
V a 



n 

(ii) For the first bound, we use (2.3.20), i. e. 

(32/3^/2)'= 



(A.i) e^/^3(^+^)/^^(y,^) = j;— ^ 



^^^r(A; + l)r((y + l-2A;)/3) 
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and the inequalities 

a + h + 2 ^ Tia + b + 2) ^ (^Tb)^ - ^' (a > 0, 6 > 0), 
which follow from a classic formula 



and 



T{a + b + 2 

af^h 



max x"'(l — x) = - 
a;e[o,i] (a + 6) 



a+b' 



1 /.1/2 /.I 

x''(l - x)'' dx > (1/2)^ / x'^ dx + (1/2)" / (1 - xf dx 



(1/2) 



Jl/2 

a+b+l ( 1 , 1 

a+1 6+1 



Break the sum in (A.l) into Si , S2 , and S3 , for {A; > 2 : (y + 1 — 2/e)/3 > 1} , 
{fe > 1 : (y + 1 - 2A;)/3 < 0} , and {fc = 0, 1 : or (y + 1 - 2fc)/3 > 1} , respectively. 
(Recah that r(0) = 00.) For Si, 

1 r(2/c/3) 1 



r(fe + l)r((y + l-2fc)/3)) r(A; + l) r((y + l-2A;)/3)r(2A;/3) 

2(2^+i)/3 r(2A;/3)r(A;/3+l) 1 



< 



r((y + l)/3) Y{k+l) r(A;/3 + l) 

2fa+i)/3 (2A;/3)2fc/3-i(A;/3)^/3 1 



r((y + l)/3) fcfc-i r(A:/3 + l) 

2^y+^)/^ (2^/3/3)^ 



-^r((y + i)/3) r(fe/3 + i)- 

Therefore 



2l.+-)/3 ^ (I^IVS)"/' ^ ,,„,3„„ 2<!«^. i,.|./o 



For S2 , we use 



r((y + i)/3^^„r(t/3 + i) """'■' ' -'nfe + ivs) " 



''<■')"<' - '-) = ihiW ^ |r(fa + i-2>.)/3)| < r(i + (2*. - . - 1)/3), 



and 



r(l + (2/c-y-l)/3) ^ 1 r(2A;/3 + 2) 



r(A; + l) - r(l + (y + l)/3) T{k + l) 

1 r(2fe/3 + 2)r(fc/3 + i) 1 



< 



r(l + (y + l)/3) r(A; + l) r(A;/3 + 1) 

1 r(fc + 3) (fc/3)'=/3(i^2fc/3)^+'^/3 

r(l + (y + l)/3) ■ r(A; + l) ■ (fe + l)'=+i 



< „./..„,. -(^.2)^'^'"/^''^ 



r(l + (y + l)/3) ^ ^ r(A:/3 + l)' 
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Therefore 

(A.3) |s,|<,(|,|»vl)^^^^2_.,HV^ 

And it is not difficult to show that 

(A.4) IS3I <b|Si| + |S2|. 

The relations (A.1)-(A.4) imply that 

D 

Proof of (2.3.23)-(2.3.24). Again we use (2.3.19). Let us choose a = ^, where 
^ = ^(y, /i) is a maximum point of 

*(a;y,/^) := -J'lna+ — + y, ae(0,oo), 

i. e. a positive root of 

(A.5) ^i'\a;y,fi)=fia + a^-^ = 0. 

a 

A root exists and is unique, since ^a(0+;y, /u) = —00, ^{oo;y, fi) = 00 and 

^l^\a;y,fx) = fx + 2a+^>0, (a > 0). 
Observe that fJ.^'^ /y is bounded away from zero. If not, then, by (A.5), 

which implies that fi^/y — > 0, contradicting y = 0{X^) = 0{fi^). 

Break the integral in (2.3.19) into /i over \t\ < /U~^'^y"'^" and I2 over \t\ > 
/i~^/^y^/^. Arguing as the part (i) of the previous proof, we bound 

\l2\<b^exp[^{C;y,fi)] J exp[-t\^ + fi/2)] dt 

(A.5) |t|>M-i/"?/i/^ 

Cexp[^{S,;y,fi)] -yV^/2 

-' (e + ^/2)V2 •« 



Turn to Ji . Since 



^^J'Hs■,y,f,) = -^ = 0{yr% 
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we have 

and 

„,3/7-l/2 

(f) 

Consequently 



/i~Cexp[*(e;y,^)] / exp f --(^ + 2^ + yt 

~^exp[^(^;y,^)] 



2^ ^ ^1^ 



,/x+2c+ye-^ 

Since y^"^ = O(^), (A.5)-(A.6) imply that A » /2, hence 

A(y,/x)~e-^'/6(2^)-i/i, 
which proves (2.3.23). D 

If we drop the condition y —> oo, then the integral in (2.3.19) is of order 



oo 

^exp[^iC;y,^i)] J exp[-t\^ + fi/2)] dt = ^exp[^(C; y,^)] 



2vr ^^/^ 



.e + ^/2 

— oo 

which proves (2.3.24). D 

Proof of (1.5). The system (1.1) is solvable iff for every cycle C of G, 

(A.7) ^ fee = 0(mod2). 

eeE{C) 

If fee £ {0, 1} are independent random variables with Pr(6e = 1) = P, the condition 
(A.7) is met with probability (1 + (1 - 2p)|C'l)/2, Kolchin [13]. 

Consider G = G{n,p = 7/n) , 7 < 1 . Let X^s denote the number of cycles of 
length s which are "bad", i. e. do not meet the condition (A.7). We need to find 
the limiting distribution of Xn = X]s>3"^"s the total number of "bad" cycles. To 
this end, observe that, with probability approaching 1 , the cycles G{n,p) may have 
are those in the unicyclic components. Let us call them u-cycles. The expected 
number of all cycles of length A; > 3 is 



,kj 2 ^ - 2k 

So 

lim lim Pr(G(n,p) has a cycle of length > A) = 0. 

A^foo ra— i-00 
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Let Yns be the total number of all u-cycles of length s. In [19] it was proven that, 
for 7 fixed, {Yns}s<A converges in distribution to {Poisson(crs)}s<^, where the 
Poissons are independent and 

£7. = , S > d. 

2s 

As 7 < 1 , we have T{'ye~^) = 7, because T{x) = xe^^-^\ for x < e~^ . Now a 
u-cycle of length s is bad with probability 

1 - (1 - 2py 

vr.s = . 



Consequently {Xns}s<A converges to {Poisson(7rcj(Ts)}s<A , whence Xn converges 
to Poisson ( J2s<A ^s<^s ) • Therefore, 



lim Prjthere are nod bad u-cycles of length A at most} = e ^s<a'^'='^^ 

n— >oo 

It remains to notice that 

1 - (1 - 2py Y 



Yl ^^^« = Yl 



2 2s 

s>3 s>3 

1 1 - 7(1 - 2p) 7. 7'.., .^ 
= — m n nil — p). 

That the same formula holds for G{n, m = "yn/2) follows then in a standard 
way. D 

Ohio State University, Columbus, Ohio, USA 

E-mail address: bgp@math.oliio-state.edu, yeumOmath.ohio-state. edu 



